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ABSTRACT 

We have found a photoevaporated disk in the Orion Nebula that includes 
a wide binary. HST/ACS observations of the proplyd 124-132 show two point- 
like sources separated by O'.'IS, or about 60 AU at the distance of Orion. The 
two sources have nearly identical / and z magnitudes. We analyze the brightest 
component. Source N, comparing the observed magnitudes with those predicted 
using a 1 Myr Baraffe/NEXTGEN isochrone with different accretion luminosities 
and extinctions. We find that a low mass (~0.04 Mq) brown dwarf ~1 Myr old 
with mass accretion rate logM ~ —10.3, typical for objects of this mass, and 
about 2 magnitudes of visual extinction provides the best fit to the data. This 
is the first observation of a circumbinary disk undergoing photoevaporation and, 
if confirmed by spectroscopic observations, the first direct detection of a wide 
substellar pair still accreting and enshrouded in its circumbinary disk. 

Subject headings: ISM: individual (Orion Nebula, M42) — planetary systems: 
protoplanetary disks — binaries: general — stars: low-mass, brown dwarfs 

^ Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. These observations are associated with program 10246. 

^present address; European Southern Observatory, Garching bei Miinchen, Germany 

■^present address: Max-Planck-Institut fiir Astronomic, Heidelberg, Germany 



- 2 - 



Introduction and Observations 



We report on the discovery of a circumbinary disk seen in silhouette against the bright 
nebular background of the Orion Nebula. Using multicolor observations taken with the 
Hubble Space Telescope as part of the HST Treasury Program on the Orion Nebula Cluster, 
we show that s ource 124-132, a ph otoevaporated disk ^1'.5 north of the Trapezium previously 
imaged in Ha (jSmith et al.ll2005l ). harbors a binary. 



Binaries represent a typical product of the gravitational collapse of cores with high angu- 
lar momentum. In fact, the majority of stars in star-forming regions are in binary or multiple 
systems (see Duchene 1999; Monin et al. 2007; and for a review Mat hieu 2000 and refer- 



ences therein). Also the ONC conta i ns a large number of binary stars (IProsser et al.lll994 



Padgett et al.lll997l : IPetr et al.lll998l : iKohler et al.ll2006l : iReipurth et all 120071 ). although the 
current counts indicate that the binary frequency, especially at the low-mass end, is lower 



by a factor 2 to 5 than in star-forming T Associations like Taurus-Auriga (IKohler et al. 



20061 : iReipurth et al.l 120071 ). Whether the relative paucity of binary systems in the ONC 



is due to the initial conditions of the cloud or to a "feedback" effect, like the ejection of 
binary companions in the cluster core caused by close dynamic encounters, is still a subject 
of debate. 

There can be three disks in a young binary system: two circumstellar disks and a cir - 
cumbinary one (ILin &: Papaloizoulll993l : lArtymowicz fc Lubowlll994l : iBate fc Bonnelllll997l ). 
with a comple x geometry due to the potential variety of alignments between the disks and 
orbital plan es (IMonin et al.l 120071). and interacting through transfer of energy and angular 
momentum (iGiinther fc Kleyll2002l ). Disks have been detected around many spectroscopic 
binaries. Direct imaging of disks around wide, well separated binary systems, either circum- 
stellar or circumbinary, is much rarer, the t wo most notab le examples being th e circumbinary 
disks around GG Tau (IDutrey et al.lll994l ) and UY Aur (IDuvert et al.lll998l ). 



The data presented here are part of the large dataset (520 images) of ACS exposures ob- 
tained for the HST Treasury Program on the Orion Nebula Cluster (GO-10246). We mapped 
with ACS /WFC an area of about 450 square arcmin nearly centered on the Trapezium Clus- 
ter in five filters: F435W (Johnson B, 420s); F555W (Johnson V, 385s); F685N (Ha+[N II] 
A6583, 340s); F775W (Cousins Jc, 385s); and F850LP (^-band, 385s). Due to the adopted 
dithering strategy, most of the field has been exposed two times so the total integration time 
is typically twice that listed. The drizzled ACS images have been visually ins pected and a 
master catalog of 219 circumstellar disks has been compiled (IRicci et al.ll2008l ). Figure 1 is 
extracted from this last paper, which also provides the fits files in electronic form. Further 
details on the HST Treasury Program observing strategy and data products will be given 
elsewhere (Robberto et al. 2008, in preparation). 
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Results 



Figure [T] shows ACS/WFC images of source 124-132. ISmith et al.l ( 120051 ) have already 
illustrated the Ha morphology of this object, a bright proplyd with a well resolved dark 
disk nearly perpendicular to the vertex of the ionization front. From their Ha analysis, they 
estimate a disk tilt angle of ^75°. The disk shows some point-like Ha emission on the west 
(right) side, which they interpret as either emission from a bipolar reflection nebula or a 
microjet perpendicular to the major disk axis. 

A more careful look at the Ha image shows that the point-like source is not centered 
on the dark disk axis, but shifted approximately 2 pixels (O'.'l) to the north, away from 
the ionization front. Since the brightness of the ionization front does not allow us to trace 
the full extension of the disk to the south, this asymmetry is probably real. The situation 
becomes clearer when our other images are considered. At shorter wavelengths (F435W and 
F555W), the point-like source remains visible to the west of the disk axes, suggesting that 
the silhouette disk is slightly tilted with the western face toward us. More interestingly, at 
longer wavelengths (F775W and F850LP), two point-like sources appear well resolved along 
the disk axes. 

In Figure [2] we plot the counts measured through the central part of the disk (column 50 
of the images in Fig. [T]). In order to remove the contribution of the sky background and of 
the proplyd we have subtracted from each cut the counts measured two pixels to the west 
(column 52). Figure [2] clearly shows the two well resolved point sources appearing at the 
same position in the F775W and F850LP filters. The F658N cut reveals the extent of the 
disk and shows that the disk is darker and possibly wider at the location of Source S. For the 
F658N filter we have also a dded a line, c orresp onding to column 51, to show the location of 



the Ha source discussed by lSmith et al.l (120051 ). Source N remains visible in the F555W and 



F435W filters, whereas Source S is perhaps barely detectable only in the F555W filter. In 
Figure [2] we also indicate the limits of the disk, showing that the disk extends further on the 
northern side. This may be due to the presence of the ionized rim, which prevents tracking 
the disk in the vicinity of the ionization front. In any case, we estimate a disk diameter of 
approximately 10 pixels, i.e., 0'.'5 (~ 200 AU), whereas the projected separation of the two 
stars is ~3 pixels, corresponding to 60 AU. 

In Table [T] we present the magnitudes of the two sources measured by integrating the 
counts in a 3 pixel diameter aperture; the sky has been measured at a close offset position 
free from ionized emission. W e have applied t he AC S zero points (Vega photometric system) 



and aperture corrections from ISirianni et al.l (120051 ). averaging the encircled energy reported 



in Table 3 of Sirianni et al. for aperture radii of 1 and 2 pixels. Being H ST diffraction 



limited, the PSF in the widest filters depends on the source color. Following ISirianni et al. 
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( 120051 ). we adjust the aperture correction in the F850LP filter to account for the extremely 
red colors of our sources. By repeating the measures with subpixel shifts of the centroid, 
we estimate an uncertainty of approximately 0.1 mag in the F775W and F850LP passbands, 
and 0.2 mag in the bluer filters due to the faintness of the source and the highly non-uniform 
background. The most striking result is the similarity between the magnitudes and I — z 
colors of the two sources, with both of them having F775W— F850LP = 1.5. 



3. Discussion 



Our images reveal the presence of a dark silhouette disk harboring two faint, similar 
objects. In the F775W and F850LP filters the sources appear clearly point-like, well sepa- 
rated, and without any evidence of bipolar extended emission from eit her sides of the disk , 
as normally observed when edge-on disks obscure their central stars (ILuhman et al.l 120071 . 
and references therein). We will thus assume that we are seeing the reddened central stars 
rather than their diffuse light scattered by the disk surface. 

To constrain the nature of the two sources, we plot their position in a color-magnitude 
diagram (Fig. [3]) and (for Sou rce N only) in a colo r-color diagram (Fig. H]) with respect to 



the 1 Myr PMS isochrone of iBaraffe et al. 



199 



calculated in o ur passbands using the 
NextGenII synthetic atmosphere models (iHauschildt et al.lll999l ). We show the stellar 



population of the ONC, we also show the displacement from the isochrone caused by extinc- 
tion and mass accretion, both combined with the synthetic stellar spectra before convolution 
with the filter passbands. For the extinction we use both the standard Ry = 3.1 reddening 
l aw and the Ry = 5.5 reddening law, co nsidered more appropriate for the Orion Nebula 
teaade fc Minkowskilll93j;_ I Johnson! 119671 ). The accretion spectrum is reproduced following 



Calvet fc GuUbrineJ Jl998h recipe, combining 3/4 of a Tes = 7000 K black body, ac- 



the 

counting for the optically thick emission from the photospheric region heated by the shock 
of accreting material, with 1/4 of optically thin emission from a dense slab (n = 10^ cm^'^) 
of photoionized gas also at Teg ~ 7000 K, accounting for the pre-shock region. The accretion 
luminosity, Laccr, and the stellar bolometric luminosity, Lphot, are both calculated integrating 
the model spectra, with Laccr finally multiplied by a normalizing factor. This approximate 
treatment is adequate for our purposes, as we are interested in the intensity in broad-band 
filters rather than in the details of the true accretion spectrum. 

Figures [3] and H] show that the observed magnitudes of Source N can be reproduced by 
different combinations of Ay, log(Laccr/ -^phot) and stellar mass. To find the best combinations 
of parameters, we calculate a large 3-d grid of models covering the full range of viable Ay, 
log(Laccr/-^phot) and stellar mass, derive the synthetic photometry and search for the best 
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match with the data. The results are summarized in Table [2] and marked as stars in Figures [3] 
and m It turns out that for each reddening law, there are typically two intersection points. 
For each reddening law, the two intersections provide combinations of lower and higher mass, 
extinction and accretion, (the low and high solutions). All four solutions place Source N in 
the brown dwarf regime. To test the robustness of this result and to find the optimal 
solution, we have assumed that the uncertainties are dominated by photometric errors and 
performed a Monte Carlo simulation. We generated a set of 10,000 trial stars populating 
a gaussian error distribution for the three magnitudes; we then iterated our procedure to 
find the corresponding mass, accretion luminosity and extinction. Figure [5] shows the results 
for the two choices of Ry- The two areas at the top and bottom of each figure correspond 
to the high and low solutions, respectively. The height of the peaks shows that the high 
solution is more compatible with Ry = 3.1, vice versa if Ry = 5.5. Only if Ry = 5.5 
the strongly degenerate high solution may lead to masses above the hydrogen burning limit 
(M ~ 0.08 Mq). 

Another clue to the nature of the sources may come from a comparison of our estimated 
mass accretion rate w ith the values normally found for similar stars. Using Eq. (2) of 



MuzeroUe et al.l (120031 ) with the stellar parameters from our Baraffe/NEXT GEN models, we 



derive for M the values reported in Tab l e El A comparison with Figure 8 of iMuzeroUe et al. 



(120031 ) ■ or Figure 2 of iMohanty et al.l (120051 ) . shows that only our low solutions provide 



mass accretion rates similar to those found in low-mass BDs, whereas the high solutions give 
values in excess by 2 orders of magnitudes. If we exclude the high values on this basis, the 
most consistent solution for Source N is the low case, corresponding to a M ~ 0.04 Mq 
brown-dwarf accreting at a rather standard M ~ 5 x 10^^^ Mq yr~^, with about Ay = 2 and 
Ry = 5.5, indicative of grains larger than those typically found in the interstellar medium. 
Similar conclusions should hold for Source S, given the nearly identical / and z magnitudes. 
Even if IR spectroscopy will be needed to firmly establish the nature of these objects, the 
possibility that the two sources have sub-stellar mass appears quite robust. If confirmed, 
this would be the first direct observation of a young brown-dwarf binary still accreting in its 
circumbinary disk. 

Whereas our images shows the presence of a large-scale circumbinary disk, the evidence 
for mass accretion points to to the presence of at least one circumstellar disk. Theory 
provides a consistent scenario, predicting that circumstellar and circumbinary disks emerge 
as the mature o utcome of the evolutio n of the original circumstellar disk in which the binary 



system formed (iGiinther fc Kleyll2002l ). Due to the exchange of angular momentum from the 
binary to the disk, an inner gap develops between the inner and outer disk regions. Material 
can flow through the gap along spiral arms, feeding the circumstellar disks and therefore 
sustaining mass accretion into the central stars. A detailed comparison of our system with 
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the extant theoretical models is hindered, however, by the fact that our circumbinary disk is 
being photoevaporated, the first time such a phenomenon is observed. The influence of UV 
flu x from an external sou rce on a circumstellar disk has been studied only for single stars 
by iRobberto et al.l (120021 ) . who have shown that the extra heating of a disk face produces 
an increase of the disk flaring angle, leading to photoevaporation and possibly a warping of 
the system. This may explain why Source S remains undetected in our bluest filters: the 
southern side of the disk is the one more directly exposed to the UV flux and the stronger 
photoevaporation may locally increase the optical depth on this side. Assuming the sources 
are identical, a line-of-sight difference AAy — 0.3™ would bring Source S below our detection 
limit at the shortest wavelengths while being compatible with the magnitudes observed in 
the red filters. 



The work of LR and ND at STScI was done under the auspices of the STScI Sum- 
mer Student Program. The authors wish to thank the anonymous referee for prompt and 
extremely valuable comments. 

Facilities: HST (ACS/WFC). 
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Table 1. Source Photometry 



Filter 


Source N 


Source S 


F435W 


23.3 




F555W 


23.0 


23.5 


F658N 


19.2 




F775W 


19.7 


19.8 


F850LP 


18.2 


18.3 



Note. — Magnitudes in Vega 
photometric system. 



Table 2. Extinction, mass, accretion luminosity and mass accretion rate for Source N. 



case Rv Ay M log(Laccr/-f'phot) logM 

(Mo) (Me yr-i) 



low 3.1 1.66 0.033 -1.76 -10.17 

5.5 2.08 0.037 -1.95 -10.26 

high 3.1 6.14 0.056 0.15 -7.83 

5.5 8.96 0.070 0.80 -7.02 
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Fig. 1.— ACS/WFC images of 124-132; Each frame shows 100 x 100 ACS/WFC pixels, 
corresponding to ~5" x 5", or ~2000 x 2000 AU at the distance of the Orion Nebula, here 



assumed to be 414 pc (jMenten et al.l 120071 ). The color picture at the bottom-right is a 
composite of the five filters with colors assigned in this way: the intensity of blue is given by 
the average between the fluxes measured in the F435W and F555W bands, the intensity of 
red by the average in the F775W and F850LP bands, and the intensity of green by the flux 
measured in the F658N filter only. Images are oriented with north up and east to the left. 
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Fig. 2. — North-south cuts (north is to the right) through the diameter of the silhouette 
disks for the 5 ACS filters (solid lines). 1 pixel (50 mas) corresponds to 21 AU at a distance 
of 420 pc. An offset has been added to the counts in each filter for illustrative purposes, 
with the horizontal dotted fines corresponding to the offset values. The horizontal dash-dot 
fine for the F658N filter is the adjacent column to the west. The position of the sources is 
marked by the vertical dashed lines, whereas the limits of the disk are indicated with vertical 
dash-dot fines. 
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Fig. 3. — Color -magnitude di a gram for the 2 point-like sources (filled circles). The solid curve 
represents the iBaraffe et al.l (Il998l ) pre-main-sequence 1 Myr isochrone, with the relative 
masses labeled along the curve. Blue dashed lines represent the displacement, for a star of a 
given mass, caused by mass accretion for increasing log(Laccr/-^phot) ratio between accretion 
and photospheric luminosity. The red dotted lines join the points of equal log(Laccr/-^phot) 
values, labeled under the 0.02 dashed curve. Dereddening vectors corresponding to 
extinction values Ay = 1, 2, 7 are shown for two choices of the reddening parameter Ry. 
The stars represent the solutions for the position of Source N compatible with the color-color 
diagram shown in Figure H] (see text). 
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Fig. 4. — Color-color diagram for Source N. Symbols are the same used in Figure 
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Fig. 5. — Distribution of mass and accretion luminosity for Source N, derived from our Monte 
Carlo simulation. The stars refer to the nominal position of the source (see text). The left 
and right panel are relative to the Ry = 3.1 and Ry = 5.5 reddening law, respectively. 



